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The FRW model with non-minimally coupled massive scalar field has been investigated in LQC 
framework. Considered form of the potential and coupling allows applications to Higgs driven 
inflation. The resulting dynamics qualitatively modifies the standard bounce paradigm in LQC 
in two ways: (i) the bounce point is no longer marked by critical matter energy density, (ii) the 
Planck scale physics features the "mexican hat" trajectory with two consecutive bounces and rapid 
expansion and recoUapse between them. Furthermore, for physically viable coupling strength and 
initial data the subsequent infiation exceeds 60 e-foldings. 
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The inflation paradigm is one of the most successful ideas allowing to explain recent precise cosmological obser- 
vations. However one of its main problems is the construction of the physically viable scenario featuring sufficiently 
long inflation epoch (> 60 c- foldings) with sufficiently large probability. In this context a lot of hope is attached to 
the models featuringa non- minimal coupling of gravity and a scalar field (e.g. the Higgs field [H-O), vector fields 01 
as well as fermions jS* (more precisely, the scalar degrees of freedom generated by the latter). Such mechanism of 
driving the inflation is also efficient in generating the correct primordial curvature perturbations. 

Among the considered scenarios the most popular one corresponds to a coupling of the form ^^ip^R, where R is 
I ■ the Ricci scalar and is a scalar field. The non-minimally coupled scalar inflaton drives the slow-roll inflation as 
^j-^ long as the non-minimal(ity of the) coupling is strong, i.e. for ^(jp ^ 1. Values of self-coupling constant of scalar 
fleld are determined by the normalization of initial inhomogeneities. For example, for the scalar field potential of the 
form V{4') = J (0^ — 2to^/A) the agreement with the observations (best flt) sets the self-coupling to ^ ^ 47000-\/A. 
In this case, the mass term does not have a signiflcant contribution either to the effective potential of the fleld, or to 
(y-^ ' the normalization of inhomogeneities. Thus, m may be of order of electroweak scale, as in the case of Higgs fleld or 
. a scalar originating from some extension of the Standard Model. Thus, one is able to explain the whole cosmological 
CO ' data by a minimal amount of new physics. Note that besides some differences pointed out in inflation from a 
\l , non-minimally coupled scalar fleld is similar to inflation in f{R) = R-\- eR^ gravity. 

— ■ At this point it is worth noting that the type of the potential considered, while usually associated with the models 
I of Higgs inflation, is not restricted just to this particular field. In fact, our studies can be easily (generalized and) 
04 ■ applied to the analysis of the inflation driven by any non-minimally coupled scalar fleld with realistic values of and 
A. 

The analysis of the running of the coupling constant A limits the allowed range of m. In case of Higgs fleld 
me (126Gey, 194GeV) with the theory error of order of 2 GeV. Thus, recent results from CMS and ATLAS 
rN . 3 {mH — 125 GeV) are both consistent with SM Higgs inflation. 

■ While the considered model is very successful on the classical level it still suffers the standard problems related with 
the presence of initial singularity, which are expected to be solved by quantum gravity. One of the leading approaches to 
provide quantum description of spacetime itself is Loop Quantum Gravity (LQG) |9l4lll| . The cosmological application 
of its symmetry reduced version, known as Loop Quantum Cosmology (LQC) *12!|, has indeed provided a qualitatively 
new picture of early Universe dynamics. The prediction of the so-called big bounce phenomenon [ISJ offered a new 
mechanism of resolving long standing cosmological problems. For example, the existence of a pre-bounce epoch of 
the Universe evolution provides an easy solution to the horizon problem, while preliminary studies indicate that the 
dynamics in the near-bounce superinflation epoch prevents the catastrophic entropy increase [l4i . J5i] usually considered 
a danger to bouncing cosmological models (following the consideration of [HI). What's even more important, the 
spacetime discreteness effects amount to a dramatic increase of the probability of inflation in the models with standard 
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potential scalar flelds [17] (see also [18|). Indeed for such models the probability of inflation with enough 



e- 



foldings to ensure consistence with 7 years WMAP data happens with probability greater than 0.999997. These results 
make the loop approach very attractive in inflationary cosmology. 
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The issue of LQC corrections to models non-minimally coupled to gravity has been analyzed in [19|, where the 
authors have introduced the LQC correction for small values of non-minimal coupling, which were not consistent with 
observations. In this paper we consider the model of FRW flat universe with scalar field coupled non-minimally to 
gravity with the realistic values of the constants A, In consequence the model features all the regimes of: strong 
{\(,\(f>^ ^ 1), medium (^^0^ > < l)j ^nd weak (|^|(/)^ ^ 1) non-minimal coupling. 

In the literature the non-minimally coupled scalar field is usually investigated via two approaches: direct calculation 
in the frame defined by the action (called Jordan frame, see e .g.l 20l') and using certain conformal transformation to 
the frame in which the field is minimally coupled to gravity [2l|. It is claimed that these two approaches lead to 
nonequivalent results [1^ [12] . In this paper we explore the idea, that the observed physics is native to Jordan frame 
(23j , however we take the point of view that it is just an emergent theory, following from the underlying "fundamental" 
theory described by Einstein frame (see sect. IIIII for more extensive justification). It is thus the latter one which is 
subject to the loop quantization in our work. 

The structure of this paper is the following: In section |lT] we introduce the model and determine the resulting equa- 
tions of motion for a scalar field non-minimally coupled to gravity. In section lTTll we perform a canonical transformation 
(to the so-called Einstein frame), which allows us to recast the theory as the one with minimally coupled matter. In 
section HVl the canonical quantization of Ashtekar-Barbero variables is carried out in the Einstein frame via methods 
of LQC giving precisely defined quantum framework. This framework is next used in section |V] as the background for 
constructing the effective semiclassical description of the system's dynamics. The resulting semiclassical equations of 
motion are then applied in section IVTl in the systematic analysis of the dynamical evolution of the universes described 
by the model. Its results are presented in section IVlIl Their general discussion involving in particular the treatment's 
limitations is provided in the concluding section rVIIII 



II. NON-MINIMALLY COUPLED SCALAR FIELD 



We start with the general (non-symmetry reduced) system of gravity non-minimally coupled to the massive scalar 
field, as defined by the following action 
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(2.1) 



where the chosen metric tensor signature is (+, — , — , — ), R is the Ricci scalar and is a scalar field (which can be 
the inflaton, Higgs field, modulus field etc). Notice that the limit of minimally coupled scalar field corresponds to 
[/(0) = 1/2. For the further studies we choose one of the most popular forms of coupling, namely 



(2.2) 



We will assume that such a non-minimally coupled scalar field (j) is the dominant matter content in the Universe. 

We next reduce the above action to the case of flat FRW spacetime g = N^At^ - a^{t)°q^ Then, after integration 
by parts the U a?'d term, one gets the reduced action 



S = 
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(2.3) 



where C/'(0) := [a^f/] ((/)). 

By varying this action with respect to a and (/>, we obtain the following equations of motion: 



iH^ + V' [- + H 



2U'd> + 2U 
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(2.4a) 
(2.4b) 



where P = 0^/2 — y is the pressure of the scalar field. In the equations above the lapse function has been set to iV = 1, 
although at this level it can safely be left unfixed. The variation 6S/SN — over it produces the (symmetry-reduced) 



^ The issue of a scalar field non-minimally coupled to gravity with non-zero spatial curvature was considered e.g. in [24| . 
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scalar constraint H = 0, where H is the Hamiltonian of the scalar field, curvature and their coupling. To define it, 
let us introduce canonical momenta of the scalar field and of the scale factor: 

dC ■ dC 

TTrf, = ^ = a^^ - Waa^, — = -UUaa - Gt/^a^- (2.5) 

d(j) oa 

For = 1 the Hamiltonian H takes the form 

H = 7T4,<p + TTad-C = -GUa^a - GJJ'^aa^ + a^SirGp, (2.6) 

where p is the energy density of the scalar field. Then the constraint H = implies 

UH^ + HU'^^^p . (2.7) 
6 

Finally, combining the eq. (|2.4all2.7|) one obtains the 2nd order equation of motion for (j> in which the curvature is felt 
only through the cosmic friction term 

. 2U'V -UV -U'^^(3U" + ^) , , 

d> + m<h = V — , 2.8 

where, as before ""' denotes the derivative over 0. 

Inflation from strong value of non-minimal coupling 

Let us introduce the slow-roll parameters defined as follows 

HH 

The universe is considered to be in the slow-roll inflation epoch whenever e, |?7| <C 1 for an extended time. For 
the coupling as in (|2.2p and the potential V{(j)) = j<j)^ these parameters can be estimated in the range of strong 
non-minimal coupling to be 

ecxr;(x^. (2.10) 

It follows then that one obtains slow- roll evolution whenever > 1, so inflation ends only when the fleld leaves the 
strong coupling regime. 

It is worth noting, that (in the considered scenario) during the inflation, <j) ^ 0(^^^/^) ^ Mpi. This implies, 
that there is no need to consider non-renormalazible terms in the potential, as here they shall be suppressed by the 
Planck scale. This is a theoretical advantage of Higgs inflation over chaotic inflation models (in which > Mpi during 
inflation). Also, despite the Standard Model Higgs fleld being a complex one, the phase of a fleld may be expressed 
by a massless scalar degree of freedom (which does not have a significant influence on the evolution of space-time) . It 
follows then that there is no qualitative difference between the dynamics of our simple model and of a more realistic 
Higgs-driven inflation. For details, see [23|. In consequence, our considerations can be applied to the physically 
relevant case 

III. TRANSFORMATION TO EINSTEIN FRAME 

An alternative approach to describing the system deflned in the previous section is its transformation to the so- 
called "Einstein frame" , where, unlike in the original physical formulation (denoted in the literature as the "Jordan 
frame" [2]), the matter is coupled minimally to gravity. This is achieved by the following (conformal) transformation 
of the metric tensor and the scalar fleld (5^,1/, (/)) {g^iv, 4>) 
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Under this change S\g^^, <j)\ = S [gfj.v,4>\ and U((j)) = 1/2, thus cj) is indeed minimally coupled to gravity (as required). 
The exphcit relation between (p and (j) may be found for all ranges of non-minimal coupling. The metric g^^ can be 
easily expressed in the FRW form in terms of t and S, defined by 

a = VWa , di = V2UAt . (3.2) 

Effective potential and energy density of transform as follows 

im^ 1 . 

The scalar constraint has the standard form of GR minimally coupled to massive scalar field, thus one can directly 
apply the known procedure of canonical quantization in the LQC framework. Ashtekar-Barbero variables in Einstein 
frame are of the form of 

sgn(u u=- 7=-5- , 6=-7VA--^, 3.4) 

where A is the so-called area gap [25| which equals A — 47r7\/3£pj [l^ and 7 is the Barbero-Immirzi parameter [26| 
equal to 7 = 0.2375 . . . ;27]. The Poisson bracket between these variables is {w, 5} = 

Before proceeding with the quantization let us note that on the classical level we have at our disposal two sets of 
variables: (u,6,0, tt^) corresponding to the Jordan frame and (w, 6, 0, tt^) in the Einstein one. Both sets are related 
with each other by canonical transformation, thus classically at the level of specified models the selection of either set 
is simply a matter of choice. This may no longer be true, once the perturbations or perturbative quantum corrections 
are introduced (20l . . On the other hand the discrete nature of loop quantization implies that its applications to 
both frames will most likely lead to inequivalent results. The issue of which choice of the frame (in general) should 
be considered the correct, physical one is under investigation p^ . although at present it remains an open question. 

Usually the selection of the set (the frame) corresponding to the physical, measured quantities is the most natural. 
In our case that would be the former set. Here however we would like to explore another (somewhat hybrid) approach: 
to consider the Jordan frame as the one corresponding to observed reality and in which all the measurable quntities 
should be calculated. However we will think about it as the "emergent" formulation, which should arise from the 
underlying theory following from the quantization in the Einstein frame. 

Such approach is motivated by two observations: one of practical, and the other of the philosophical (or rather 
aesthetical) nature. First, the Einstein frame is a much more natural candidate for quantization. Indeed, in that frame 
gravity and matter fields are minimally coupled and the gravitational part of the action takes the standard GR form 
(Einstein-Hilbert action term). Thus, the quantization procedure is natural and well understood. Second, the classical 
structure (variables) underlying the loop quantization originates from classical GR and encodes its symmetries. Thus, 
provided that loop quantization is applicable to the considered system at all, Einstein frame is the natural arena for 
it. 

We employ the above idea via implementing the following procedure for the quantization of the system and the 
description of its dynamics: 

(i) First, we provide the complete quantum description in Einstein frame. 

(ii) Next, the quantum framework is used as the basis for constructing the effective description, used in turn to 
analyze the dynamical behavior. 

(iii) Finally, the results of dynamics are translated to the "physical" Jordan frame. 



IV. LOOP QUANTIZATION 

As discussed above, we quantize the system in Einstein frame, choosing the polymer representation for the geometry 
degrees of freedom and the Schrodinger one for the matter ones. Thus, the quantization pro cedure is a full analogy 
of the one applied to the FRW cosmology with minimally coupled massless scalar field in [25[ . 



^ We use variable b with opposite sign than usual in order to regard v as our configuration variable. 
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To start with, we employ (the initial part of) the Dirac program, first quantizing the system on the kinematical 
level (ignoring the constraint). The resulting kinematical Hilbert space, Hkim is of the form 

where M is the Bohr compactification of the real line. A convenient basis of Hgi is formed by the eigenstates of the 

~ C-3 

oriented volume operator V = d 

V\v) = 27r7\/A4iw|w) =: a\v) , (4.2) 

where \v) is eigenstate of operator v with eigenvalue v. They are orthonormal with respect to Kronecker delta, that 

is {v\v') = 6v,it'- On the other hand, the basis of He/, is provided by the generalized eigenstates (0| of the operator 
- a quantum counterpart of the field value 0: 

VxeH^ (^|0-^/|x) =0. (4.3) 

At this point a short comment regarding the exact physical meaning of the quantity V is in order. In order to obtain 
a meaningful canonical description of an isotropic model with noncompact spatial slices (which is the case here), one 
introduces into the theory an infrared regulator - some fiducial cell V constant in comoving coordinates. The quantity 

V — d^ here corresponds exactly to the (oriented) volume of that cell. Of course, one has to remember that in order 
to obtain consistent description one has to make sure that the theory admits a well defined regulator-removal limit. 

Given "Hkin, we now select the basic operators defined on some dense domain of it. Two of them are the operators 

V and (j) (|4.21 14. 3|) . The remaining two are the unit shift operator N and the scalar field momentum n^: 

N\v) = \v + l}, ^i = ^^d^- (4.4) 

The pair w, N is the equivalent of the operators of quantum flux and holonomy in full LQG, where the holonomy-fiux 
algebra provides the basis for quantization 0. particular, i represents (the power of) the fiux across the unit 

surface (in comoving coordinates), whereas N the holonomy along straight line [281] . 

The set (-5, iV, </>, tt^) is suflacient to construct the quantum counterpart of the Hamiltonian constraint - the next 
step in Dirac program. The details of its construction are presented in [2^1. However, here we choose a bit different 
and more convenient symmetric factor ordering, the one used in [29| . The resulting operator is 

Q_-/^ I - ' I 

H = Hgr®l0 + H0, %, = —J\i\{N^-N-^fJ\il (4.5a) 

o<y. 

H0 = ^^i7r? + ||5|n0), (4.5b) 

where a is defined via (|4.2I) and V{4)) is the effective potential given by p.3p . 

Having at our disposal both "Hkin and the quantum constraint H we can complete the Dirac program, constructing 
the physical Hilbert space (out of states annihilated by H) via group-averaging procedure [30| , and introducing the 
family of unitarily related observables to capture nontrivial information about the evolution of the described system. 
Such observables can be again constructed out of self-adjoint operators on Hkin (kinematical observables) via group 
averaging. In generic situations, however, the form of the resulting physical observables is very different than the 
original kinematical ones. Therefore, to obtain a technically manageable theory, one usually resorts to the so-called 
deparametrization, introducing into the system a suitable matter field and using it as an internal clock. This is exactly 
the method originally used in psl ]. 

One such deparametrization, particularly successful in full LQG, has been introduced in [sij. There, the time 
variable is provided by irrotational dust (see [slj for studies of dust frames and [s^ for application of full dust frame 
to LQG). Then, the synthesis of three elements (the specific matter field, the fixing of the time gauge provided by 
the proper time of dust "particles" , and the diffeomorphism-invariant formalism of standard LQG) allow to build 
a technically manageable completion of gravity quantization program. Its reduction to isotropic cosmology has 
been presented in [23]. Its main properties are: (i) the Hamiltonian constraint becomes true Hamiltonian, (ii) the 
kinematical Hilbert space Hkin becomes (precisely) the physical one, so do the kinematical observables, and (Hi) the 
evolution of the physical state is described by time-dependent Schrodinger equation 



-ihdf'i'iv,^) =H*(S,(^) , 



(4.6) 
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where the time i provided by the gauge choice agrees with the cosmic time and the Hamiltonian H (in our case 
given by (j4.5p ) is self-adjoint for any system of gravity coupled to non-exotic matter [34| . This quantum description 
provides an excellent basis for further analysis of the dynamics. It can be performed on the genuine quantum level, 
although in this paper we will employ the so called effective description which will be introduced in the next section. 
In order to describe the dynamics we select the following set of obscrvables, convenient in cosmology: scaled 

quantum volume |{;|, operator sin(5) :— (i/2)[iV^ — -^^^], field (f) and its momentum ttt. All these operators, being 
originally kinematical observables in irrotational dust deparametrization automatically become physical. 

Another observable especially relevant for cosmology and in particular for our further studies is an operator corre- 
sponding to Hubble parameter. It takes the form: 

H^~-^[N^^N-% (4.7) 

At this point it is worth noting that the dust field, although convenient, is not the only possible choice of the 
internal time. One alternative, viable in the cosmological context is the choice of additional - massless scalar field. 
Such field has been used exactly in the pioneering work [2'5| and is in fact the m ost p opular choice of clock in LQC. 
Its application to the full theory is well defined (as variation of the formalism of [35|) at least formally, however the 
complications related with the properties of the scalar field keep its application to non-cosmological settings beyond 
current technical reach. In the homogeneous cosmology context it is manageable, although not free from (minor) 
technical difficulties, such as: (i) evolution is provided by Klein-Gordon equation rather than Schrodinger one, (m) 
in presence of any non-negative scalar field potential the Hamiltonian loses self-adjointness, and {Hi) the physical 
Hilbert space is a proper subspace of "Hkin and its identification requires explicit knowledge of spectral properties of 
H. We discuss this formalism in more detail in appendix \K\ 



V. EFFECTIVE THEORY 



At this point we have at our disposal the physical Hilbert space, time, true Hamiltonian operator and the set of 
physical quantum observables. These components are sufficient to systematically study the dynamics of our system. 
The analysis on the genuine quantum level can be performed by the methods of [36j . However, since the aim of this 
paper is to obtain maximally complete picture of the semiclassical sector of the theory we resort to effective methods, 
leaving the verification of the results against the genuine quantum ones to future work. To arrive to such effective 
description, we employ the method introduced in |37| (allowing in principle to account for arbitrary order quantum 
corrections, see e.g. 38]) in its 0th order. In this order the resulting set of equations of motion is equivalent to the 
one provided by heuristic methods of the so called effective dynamics of LQC [39J. Its comparison against the genuine 
quantum dynamics in many LQC systems [1^ l40l - |42j has shown that it mimics the full quantum evolution to high 
level of accuracy. Furthermore, in some cases this result has been confirmed analytically [43!| . 

Technically, the method reduces to replacing the basic operators in quantum Hamiltonian by their expectation 
values evaluated on the semiclassical (sharply peaked) statefo. In our case, this is 

^"h^(^)", TVt^ (iV)e'^/^ (5.1) 
Applying this mapping to (|4.5p . we get the following "classical" effective Hamiltonian 

2 

= -^|e|sin2(&) + ^ + ^y{i) = Eos. (5.2) 

where Ecs is the energy of the dust clock field. The full set of the effective equations of motion is then provided by 
Hamilton's equations 

2 

dv 6nG , , db 3nG '^x 2a ~ , , 

|w| sin(6) cos(6) -= = — - sin2(5) + — ^ - -^V (5.3a) 



di ah di ah ahv^ ft^ 

d(f) _ TT^ dTT^ ^ a\v\ dV 

di a\v\ di h di 



(5.3b) 



There we assume implicitly, that there exists a sufficiently large space of such states. Verification of this assumption requires testing 
the dynamics on the genuine quantum level. 
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In order to eliminate the effects of the dust field on the dynamics, we set 

= 0, (5.4) 

thus reducing the presence of dust to a "dust vacuum" . Under this condition the set of equations of motion ()5.3a|) 
and (|5.2l 15. 4|) becomes equivalent to the following set 



3i/' = SvrGp 1 - — , — - = -47rG p + P 1-2-^ , (5.5) 

\ Per J dt V / \ p^rj 

where dv /dt = SvH and pcr = S/SttGAj'^. The Big Bounce in Einstein frame appears for p = p^r- 
It is worth noting that (j5.5l) can be converted to Jordan frame via (j3.ll) , giving 

6f/i^^ + 6^'i^^ + ^^f^.8.G(. + ^^f!) + , (5.6a) 



Note that for S^cjP ^ 1, unlike in the case of LQC with minimal coupling, at the moment of the Bounce we have 

p > Per- 

The equation of motion for cj) is of the form (|2.4ap : surprisingly, the effective rolling force of the scalar field is not 
changed by the LQC correction. Thus, the LQC correction does not influence the classical equation of motion for 0, 
which (after the substitution of = gives eq. 



VI. THE ANALYSIS OF THE DYNAMICS 

The set of equations of motion in the Einstein frame (|5.3p allows in principle to evaluate the time evolution of 
the canonical data {v,b,(j),n^) up to the caveat that the form of (j){(l>) is needed to specify V{(j)). This dependence 
is provided through the differential relation p.ip . Therefore, for the sake of precision, it is much more convenient 
to formulate and integrate the mixed set of equations of motion for the variables (v, b, cj), tt^), where the equation for 
d(t>/dt follows from (j5.3bp and (|3.1|) . This set allows to determine all the relevant physical parameters. For technical 
reasons (faster evolution), it is more convenient to evolve the system with respect to the "physical" time t of the 
Jordan frame. Using (|5.3|) . we arrive to the final set of the evolution equations: 



dv 6ttG I ; — -, , db 



V^U{^)\v\ sin(6) cos(&) — = ^2Ui(t,) 



dt ah dt 



^sm^{b) + ^^-—V 
ah ahv^ h^ 



4[/3((?!)) TT^ d-K^ _ a\v\ / 4[/3(0) dV 



dt V C/(0) + 3f/'2(,/)) a\v\ dt h V [/((/)) + 3[/'2(,/)) dfj) ' 



(6.1a) 



(6.1b) 



In order to specify the initial conditions for the time evolution, we exploit the fact that each physical trajectory has 
a distinguished point, namely the big bounce in the Einstein frame (at which b = 7i'/2). We chose it as our "initial 
point" , setting t = there and evolving the initial data (set there) both forward and backward in time. Since the set 
(j6.ip is homogeneous in v, |f| we have the freedom in choosing its initial value. We set it to 1. Having the geometry 
degrees of freedom set, we fix the scalar field momentum tt^ as the function 7T^{v,b, cj), E^g) using the Hamiltonian 
(|5.2p . To eliminate the influence of the dust "clock" field on the dynamics we further set its energy EcS to zero. As 
a consequence we have at our disposal a 1-parameter family of initial data (t — 0,v = l,b ^ 7r/2,0i„,7r^ = n^{v = 

l,b — 7r/2, (pin, Ecs — 0)) parametrized by the value <^in of the scalar field at the bounce point in the Einstein frame. 

These initial data, supplied (and determined) by the set of constants ^, A, was then evolved using the fifth order 
adaptive Runge-Kutta (Cash-Carp) method using the effective dynamics module of the Numerical LQC library de- 
veloped by T. Pawlowski and J. Olmedo. The raw results of the simulations were next postprocessed with the use of 
Mathematica and gnuplot software. 



* One needs to remember that v encodes the information about the physical volume of the chosen region of spacetime, whose choice for 
noncompact spatial topology of the universe is arbitrary. 
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In actual simulations the value of A was set to 1/2, which is the same order of magnitude as in the Standard 
Model. Due to the running of the coupling constant one cannot predict the precise value of A around the Planck scale. 
However recent results from ATLAS and CMS [1| suggest that A shall not be too big, since the mass of the Higgs is 
close to the minimal allowed value which does not violate electroweak vacuum. Thus, the value of A assumed by us 
is realistic. 

The relation between coupling constants ^ and A is given by the normalization of primordial inhomogeneities, 
which gives ^ = 47000-\/A. However, to confirm the robustness of the results and further analyze the qualitative 
behavior of the system, several different values have been considered: the sequence of lower values ^„ = {4.7 x 
WVX; n = 0, 3}. To test the wide range of trajectories the initial scalar field values were chosen from the interval 
0in S (— 10^~^/^, 10^~^/^). For most simulations we selected 20 initial points distributed unformly within this interval. 
These initial data have been next evolved till the time ifi„ — 2 x 10'^, 10^ depending on the simulation (and, in 
particular, on the constant ^). The results of these simulations are discussed in the next section. 



VII. THE RESULTS 



The results of our studies can be summarized in the following set of points: 

• In the (underlying) Einstein frame we observe the standard for LQC picture of a single bounce separating 
contracting and expanding epochs of the universe evolution. The "bare" matter energy density p and Hubble 
parameter are bounded by their respective critical values pc ~ 0.41ppi and He ~ 1.3£p^. 

• In the Jordan frame, here conjectured to represent the observed dynamics, the standard bounce paradigm is 
slightly changed. The single large energy density epoch still separates two long epochs of contraction and 
expansion. The process of the bounce itself is however modified: the evolution of the scale factor features the so 
called "mexican hat" shape (see Fig.|4]) - the sequence of bounce, ultrarapid nonadiabatic expansion ending with 
recollapse, similar nonadiabatic epoch of contraction and the final (second) bounce, after which the universe 
expands to the classical regime (see also Fig. |3] for the Hubble parameter evolution). The time between the 
bounces is of the order of Planck time. 

• Outside of the "mexican hat bounce" the dynamical trajectory approaches quickly the one predicted by GR. 
In particular in the future of the bounce the value of the field </> grows to certain (depending on the trajectory) 
maximal value at which point the slow roll inflation starts. Due to the time symmetry of the equations of 
motion the inflation after the bounce is accompanied by the "slow roll" deflation before it. The initial Hubble 



parameter at the onset of inflation is proportional to (pi-nn (see Fig. 2(a) I 



The number of e- foldings during the inflation is estimated to be N+ ~ 3/16£0^ and for the physical value of 
^ ss 4.7 • 10^ exceeds 60 within all the studied range of initial data (see Fig. 2(b) I. The same estimate holds for 



the deflation: N- ~ S/lG^^^^jg, where ^dofi is the minimal value of the fleld (reached before the bounce) 

Within the precision of our estimates, the product 7V_|_7V_ of the e-foldings during inflation and deflation does 
not depend on the initial data. It is the function of ^ and A only, in particular growing with ^. For the selected 

1 2 



values of ^, A it equals N^N^ = 1 4.4 • 10 I (see Fig. 2(b) I. This also implies that to have an inflation with 



N < 60, there must be a deflation with iV > 10^. This is an extremely asymmetric situation, and therefore 
highly unlikely at least on the intuitive level: we thus expect that the inflation with > 60 is extremely 
probable - with much higher probability than in the case of the standard chaotic inflation. 

• At the late time each trajectory reaches the standard inflationary attractor and the reheating (see Fig. [5] for the 
illustration on the example of the unphysical model corresponding to ^ = AlVX). 



VIII. CONCLUSIONS 



We investigated the model of non-minimally coupled scalar field coupled to gravity via additional term in the 
Lagrangian density: ^£,(j)^R and with potential jXtj)"^- This system has been studied via effective dynamics methods 
within loop quantum cosmology framework. In the process of constructing the treatment we implemented the idea 
where the so-called Jordan frame represents the observed (physical) dynamics, however it is an emergent framework 
coming from underlying one represented by Einstein frame. The studies have shown that, similarly to other models 
within LQC, the high density region connects two semiclassical branches of the universe: contracting and expanding. 
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FIG. 1. Global view 1(a) and the zoom into particular sectors of the phase portrait with respect to matter degrees of freedom 

shows 



1 c 



for the physical case of C = 47000\/A. In [T(b)] we zoom the region around the bounce in the Einstein frame. Figure 
the trajectories at the onset of inflation (at (p = 0), providing the initial conditions for inflation. We do n ot sh ow the attractor 
trajectory, as the length of inflation exceeeds the technical (time) limits of our numerical simulations. 1(d) shows the same 
sector as |l(b)| in the logarithmic scale to visualize the behavior of trajectories for small (j>. 



However the process of the bounce itself is modified, featuring the "mexican hat" shape of the scale evolution with 
highly non-adiabatic epoch of expansion— >contraction between two bounces. The expanding epoch features the slow 
roll infiation whose length for all the investigated trajectories exceeds 60 e-foldings. These results, although solid 
within the specified model, cannot be treated as final due to the limitations of the applied model itself. We discuss 
these limitations below. 

First, the quantization procedure applied here assumes the role of Jordan frame as the observed one and the Einstein 
frame as the underlying (somehow more fundamental) one. Such approach raises a (natural for non-minimally coupled 
systems) concern, that the description may admit a Lorentz symmetry violation, which in principle could imply 
disagreement with the cosmological observations (for example, through the structure of the primordial gravitational 
waves). To the best of our knowledge there does not exist any definite result indicating that this would be the case 
but the issued in itself remains still open. 

Furthermore, the quantization procedure in itself requires that the gravitational part of the Hamiltonian can be 
expressed in the form where matter is minimally coupled to gravity. However, this assumption cannot be satisfied 
for a vector field non-minimally coupled to gravity, as the the transformation described in the eq. p.ip would still 
leave non-minimal coupling terms. This problem cannot be solved by the transformation to Einstein frame even if a 
non-minimally coupled vector field is subdominant, i.e. when the presence of a vector field does not break isotropy 
of the background metric tensor. In these cases, one needs to quantize the system directly in Jordan frame. This is 
especially relevant for calculation of the quantum gravity effects for models such as vector infiation [4|] . 

For any multi-inflationary scenarios with at least two fields non-minimally coupled to gravity, one can transform 
the action to the Einstein frame, obtaining the standard curvature term in the Lagrangian density. However, after 
this transformation, matter fields in Einstein frame would feature non-standard kinetic terms. In that situation one 
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(b) 



FIG. 2. 2(a) Projection of the phase portrait onto the . 



H plane. The trajectories reach the onset of inflation (marked by blue 
crosses) which points lie on a straight line (as expected in the strong non-minimal coupling case, where H ~ ±-,/\JV2i\(j>\). [2(b) 
shows the number of e-foldings during the slow-roll inflation after the bounces (red line) and the deflation before the bounces 
(blue line) as a function of initial value of the field (labeling the trajectories). Notice that even the trajectories with very high 
time reversal asymmetry (considered to be unlikely) give more than 60 e-foldings during inflation. 




FIG. 3. Evolution of the Hubble parameter as a function of time (in Planck units) for the symmetric trajectory 
and for two non-symmetric ones. Fig. 3(a) shows the general view around the epoch of high energy density. In 3(b) 



<j}iu = 0) 
we zoom 

around t = Q for the symmetric trajectory only (with the recoUapse at t = 0), to show that the transition is in fact smooth. 
3(c) visualizes two bounces in Jordan frame. 3(d) shows the regions of the onset of inflation and the end of deflation, where H 



is approximately constant. 



has to resort to the quantization schemes of non-canonical kinetic terms used in QFT. 
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FIG. 4. The evolution of the scale factor for the symmetric trajectory and for two non-symmetric ones 



±0.0055) |4(a 



and its close-up around t — 4(b) showing the characteristic "mexican hat" shape with two bounces and recoUapse between 
them. Note that for non-symmetric trajectories the "tip of the hat" (i.e., the moment of recoUapse) is shifted in time and less 
evident. 
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FIG. 5. Phase space (matter sector: cj>-ij> plane) and evolution of Hubble parameter for the ^ — 47\/A toy model. For such 
a small value of the coupling we are able to evolve the system through the whole inflation. Fig. 5(a) shows the global view 
of th e 4> half. 5(b) shows the approach of the trajectories towards the attractor (ending in the reh eating at <^ ~ \/A<^'^). In 
c) we present the time evolution of the Hubble parameter during the inflation/deflation epoch. 5(d) is a close-up of the 
same evolution around reheating. The qualitative features demonstrated here are expected to be also present in the realistic 
C = 47000VA case. 



12 



One has also to remember that the results presented in this paper have been obtained via the 0th order effective 
dynamics. This choice of methodology is motivated on the one hand by the reliability of the method itself proven for 
other models of LQC and, on the other hand, by the aim of the study: verifying the robustness of sufficiently long 
inflation in the model considered (which in turn required probing maximally wide region of the space of solutions). 
However, the reported results can be treated as preliminary only, since the state-dependent parameters may affect 
the dynamics especially in the strongly nonadiabatic epoch between the bounces. Therefore, the effective method 
has to be confirmed by a higher order effective treatment of [37| and especially by the studies of the dynamics on 
the genuine quantum level within the framework of (29| (the latter in particular to determine the viable cutoff of the 
effective treatment order). These studies will be performed in future work. 

Finally, let us note that the modification of the dynamics in the high energy epoch may also affect the structure of 
the inhomogeneities in the universe. In particular, the possible formation of shock waves in the strongly nonadiabatic 
epoch between bounces may in principle affect the structure formation in the expanding universe. Addressing these 
issues requires however extending the simple isotropic formalism used in this paper to inhomogeneous situations. The 
first step in this direction is the formulation of a well defined canonical formalism (in particular the adapted Ashtekar- 
Barbero variables) directly in the Jordan frame, followed next by the quantization within LQG framework. Then, 
an application of recent developments in the deparametrization techniques of LQG 31j will allow to unambiguously 
control the dynamical sector of the full theory. 

This in turn will allow to construct a precise perturbative framework. Methods of building such framework are 
already available in the context of LQC. They are: the hybrid quantization 14411 . strictly implementing the unitary 
evolution, or Quantum Field Theory on Quantum Cosmological Spacetime |45| . already showing some preliminary 
successes in determining the primordial perturbations [4^. Applying these tools to the scenario studied here will 
ultimately allow to test its consistency with CMB observations and identify possible effects characteristic for it. 

ACKNOWLEDGMENTS 

We would like to thank Mikhail Shaposhnikov and Dominika Konikowska for helpful discussions. The work was 
supported in parts by grants of Minister Nauki i Szkolnictwa Wyzszego no. 182/N-QGG/2008/0, N202 091839 and 
N202 104838. The bulk of numerical simulations has been performed with use of the Numerical LQC library currently 
developed by T. Pawlowski and J. Olmedo. 



Appendix A: Deparametrization via massless scalar field 

One of the steps performed in order to arrive to an effective description of the considered system in Sec. |V] was 
the introduction of the dust matter field and next the deparametrization of the system with respect to it. While 
this particular choice of the clock allows the full formal control over each step of the construction of the effective 
description, it is by no means the only one possible. Indeed, deparametrization can be performed (at the cost of a 
limited loss of precision) with respect to other matter fields of sufficiently good properties, for example the massless 
scalar considered in the original LQC works or the (appropriately tailored) population of Maxwell fields [13] . As 
long as we restrict the expansion to the 0th order and assume vanishing mean energy of the clock field, the formalism 
leads to exactly the same equation of motion (j5.3|) as in the case of dust clock. In this appendix we briefiy outline 
the construction of the effective description with use of the massless scalar field. The exact same procedure can be 
applied to Maxweh fields of 0|. 

Our starting point is the modification of the classical Hamiltonian constraint in the Einstein frame, namely adding 
to it the massless scalar field term: 

H ^ H' = H + , = , (Al) 

2a\v\ 

where p^p is the canonical momentum of the clock field ip and the variable v (and the constant a) is defined via 
eq. (|4.2p . This constraint is next deparametrized with respect to in effect being brought to the form 

pI = Q = -2a\v\R , (A2) 



^ See [25I of the detailed description in the case of vacuum with massless scalar field. 
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which in turn can be understood of as representing the free system evolving with respect to (f. Applying polymer 
quantization as in Sec. IIVI (see (25| for further information) we arrive to the generalized Klein-Gordon equation 

-dl^{v,^) = m{v,cj,) , e = -2a:|«|H:, (A3) 

where : • : means the symmetric ordering. Action of on a dense domain in 'Hkin reads (following (4l| for its 
gravitational part) 



0gr 



e = Ogr «) l0 - 6^ , (A4a) 



|{)|(7V-iV~i)VH , 00= Igr «) TT? + 2a2v2 (g) F(<^) . (A4b) 



Taking one of the superselected sectors (e.g. the one corresponding to the positive frequencies), one arrives to the 
formal Schrodinger equation 



id^^{v, 0) = V m{v, (/)) . (A5) 

To write it precisely one needs to overcome several obstacles as 9 is neither positive definite nor essentially self-adjoint 
jTsi] . since it contains the potential term acting as effective positive cosmological constant on = const plane. It is 
however expected to admit a family of non-unique selfadjoint extensions. Therefore one has to choose one extension 
QjS (where /3 is an abstract extension label). Then, upon fixing the gauge t — ip, the Schroedinger equation can be 
written as 



tdt^iv, 4>) = H^vl;(^;, cj,) := ^ [e^] + vI/(«, 4>) , (A6) 

where [•]+ denotes the positive part of the operator. Since only the positive part of 9^ enters the Hamiltonian, the 
physical Hilbert space does not coincide with T-Lmh- Instead, it is its proper subspace Hp spanned by the eigenfunctions 
of [9^] + . In principle, different extensions could produce different physical subspaces, as it is the case for positive 
cosmological constant [^l]. The consequence of the latter is that every operator defined on 'Hkin has to be further 
projected on Hp to define the physical operator, thus the operators like v lose their simple analytical form. 

Proceeding with the construction of the effective description requires a heuristic assumption: namely, that for 
sufficiently large class of states the projection does not affect the action of the elementary operators significantly. 
Given that, we can replace v" and N by their expectation values, arriving to the Hamiltonian 



Hcff = ^SttGv^ sm\b) - TT? - 2a2„2^(0) = =; . (A7) 

Similarly to the dust time case, this (also heuristic) step corresponds to the 0th order effective expansion introduced 
in [s^]- Here, however, this procedure is not fully justifiable, since it ignores the effects of inequivalent self-adjoint 



extensions on the dynamics (see the discussion in (38() - effects that, although small, are usually present [42|. 

If we accept the heuristic steps discussed above, we arrive to an effective model described by the Hamiltonian (|A7p . 
To compare it against the treatment of Sec. |V]we note that: 

(i) For Pip ^ the model is equivalent to the one described by the Hamiltonian 

H^g = SttGv"^ sin2(&) - TT? - 2a^v^V{4>) . (A8) 

The equations of motion of this model are equivalent to the one of (jA7p with the lapse function rescaled via 
TV ^ 2p^iV. 

(ii) Under the change of sign and the lapse function transformation N -> N/ {2av) the Hamiltonian (|A8p transforms 
into the effective Hamiltonian of the system with dust time (|5.2I) . 

In consequence, at the level of 0th order effective dynamics the deparametrization with respect to massless scalar field 
and irrotational dust are physically equivalent. 
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